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IONIZATION OF NO AT HIGH TEMPERATURE* 


C. Frederick Hansen** 


INTRODUCTION 

Space vehicles flying through the atmosphere at high speed are 
known to excite a complex set of chemical reactions in the 
atmospheric gases, ranging from simple vibrational excitation to 
dissociation, atom exchange, electronic excitation, ionization, and 
charge exchange. The major reactions in high temperature air were 
measured in the 1960 's, generally in shock tubes, at temperatures 
up to about 10000°K ; and Arrhenius type equations were then fit to 
the observed rate coefficients. Although rate data is typically 
rather scattered, even for measurements made with a fixed 
experimental arrangement — and different experiments often yield 
widely different results — yet for some engineering purposes the 
uncertainties are not excessive. This occurs because endothermic 
rate coefficients, at least, change so rapidly with temperature 
that in flow field calculations where temperature gradients are 
high, as they are in the blunt body regions or boundary layer 
regions of high velocity vehicle flow, the uncertainties in the 
rate coefficient have a minor effect on species concentration and 
thermodynamic properties as a function of distance along the flow. 

However, the exothermic rate coefficients can introduce 
significant uncertainties in cooler regions of flow such as the 
inner boundary layers adjacent to the body surface or in the 
expansion flow along the aft part of the vehicle and in the 
vehicle's wake. Moreover, the rate coefficients are now needed at 
much higher temperature and lower density than formerly , for 
application to the new class of space vehicles planned, such as the 
AOTV — Aeroassisted Orbit Transfer Vehicle. Temperatures up to 
40000°K are well beyond the range that seems accessible to 
experiment and extrapolation of the formulae which fit the lower 
temperature data is uncertain unless a theoretical model can be 
devised to justify the extrapolation. Finally, the new class of 
space vehicles will decelerate at very high altitude (more than 60 
km. above the earth's surface) where the density is so low that the 
shock heated gases will often not be in either thermal or chemical 
equilibrium. Since most of the measured reaction rates have, been 
made at conditions where equilibrium is believed to obtain, a 
theoretical model is also desired which can treat particularly the 
effects of vibrational and electronic nonequilibrium on the rate 
processes . 

* The purpose of this paper is to develop simple arguments for 

the temperature dependence of the reactions leading to ionization 
of NO, including the effect of vibrational-electronic thermal 
nonequilibrium. NO ionization is the most important source of 
electrons at intermediate temperatures and at higher temperatures 
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provides the trigger electrons that ionize atoms. Based on these 
arguments, recommendations will be made for formulae which fit 
observed experimental results, and which include a dependence on 
both a heavy particle temperature T and different vibration- 
electron temperatures T v and T e . In addition, these expressions 
will presumably provide the most reliable extrapolation of 
experimental results to much higher temperatures. 


TEMPERATURE DEPENDENCE OF EXOTHERMIC AND ENDOTHERMIC REACTIONS 

The rate coefficient for a two-body collision induced reaction 
can be expressed guite rigorously as the reaction cross section S 
times the collision velocity integrated over a Maxwell-Boltzmann 
distribution of collision energies 1 . 

oo 

k - — ( Six) x e x dx Eq.(l) 

s J 


where u is the mean velocity of the collision pair in center of 
mass coordinates, s is the symmetry number unity for dissimilar 
particles, two for similar particles, and x is the dimensionless 
collision energy in units of kT (i.e. x = E/kT = /xu 2 /2kT) . 

If the collision is exothermic with no activation energy, the 
cross section is expected to be a relatively slowly varying 
function of energy. For example, if the cross section varies 
inversely with the m th power of collision energy 

5 = _£° - (kT)~ m Eq . ( 2 ) 
E m x m 


then the integral of Eq.(l) is analytic if m < 2 

k = 1 - g ° r (1-in) « T~ {m ~ 1/2) Eq . ( 3 ) 

exo ( kT ) m 


However, when the reaction is endothermic, the integration is 
performed over x only from x 0 to infinity, or over all y = (x - x D ) , 
the dimensionless energy in excess of the threshold energy x 0 . 

k - -H 


/ siy) 


(x+y) e~ y dy 


Eq . ( 4 ) 


Eq . ( 4 ) is the form of the usual Arrhenius expression for 
endothermic reaction 
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k - d.T)e- EJn Eq ‘ (5) 


where the temperature dependence of the preexponential factor 
depends on the average collision velocity u and the definite 
integral determined by the shape of the cross section function S(y) 

OO 

C(T) - — f S{y) (x +y) e~ y dy Eq.(5a) 

S J 


An additional factor that needs be considered is that while 
highly excited electronic states do not contribute much to 
exothermic reaction because their population numbers decrease so 
rapidly, the endothermic reactions to not usually occur from the 
ground state but from excited electronic states that have 
activation energies the order of kT. The exponential decrease in 
population number of the excited state is balanced by the 
exponential increase in the rate coefficient due to the lower 
activation energy. The total rate coefficient is the sum of rate 
coefficients from each state i weighted by the probability that the 
i«* state is involved in the collision. This remains an expression 
like Eg. (5) with the same activation energy E 0 


k 



- \ T,9,C 1 iT)e 


{E 0 -E i )/kT e -EjkT 



( E^tn) 


e -E a /kT 


Eg. (5b) 


where n t and n are respectively the densities of the reactive 
molecule in the i“ state and the total number density of all 
states. The preexponential term is a sum of all the preexponential 
terms Ci(T) weighted with the degeneracies g t divided by the 
partition function Q; these sums may include all the rotational and 
vibrational states of internal energy as well as electronic states. 

The temperature dependence of endothermic reaction is 
dominated by the exponential factor exp(-E 0 /kT) whenever the 
activation energy E a is large compared with kT. Then it is not 
possible to arrive at a good estimate for the temperature 
dependence of the preexponential term from experimental data 
because the slope of an Arrhenius plot of ln(k) as a function of 
1/kT is very nearly —E,, no matter what preexponential function is 
chosen. For example, an endothermic rate coefficient that has been 
widely used 5 ' 3 for the heavy particle impact dissociation of 0 2 is 
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_ 3 . 7 4-xlO 18 e -594oo/r cc/mol-sec Eq.(5c) 

T 


However, as Fig.(l) shows, the slope of the Arrhenius plot is 
changed very little if the exponent on the preexponential 
temperature terra is changed ±1/2. This particular example is 
chosen to illustrate the point because the data provided by Camac 
and Vaughan 2 on 0 2 dissociation by Ar impact is about the least 
scattered data for high temperature endothermic reaction that has 
been obtained. Still enough scatter remains so that the least 
squares deviation is almost the same for any one of the three 
functions shown. In most cases the scatter in measured rate 
coefficients is much greater; thus the best way to arrive at the 
correct temperature dependence of rate coefficients is to determine 
the exothermic rate, if possible, and then calculate the 
endothermic rate with the equilibrium constant, which is a well 
known expression involving the partition functions of the 
chemically reacting species. 

To develop the analysis further , the shape of the reactive 
cross section as a function of collision energy must be given, 
either as an experimentally measured parameter, or as a theoretical 
estimate . 


REACTION CROSS SECTION THEORY 

Various empirical or semiempirical approximations have been 
tried for the reactive cross section. Qualitatively the cross 
section must vanish at collision velocities less than the 
activation energy, then it increases rather rapidly as collision 
energy exceeds the threshold, maximizes at something less than the 
total scattering cross section, and finally decreases at much 
higher collision energies. Hansen 4 has suggested a cross section 
function for collision induced dissociation of 0 3 , N 2 , and NO which 
increases as the square of the excess collision energy above 
threshold and eventually decreases as 1/E 


— - 2(1 -Ej E) 2 {Ej E) Eq . ( 6 ) 

S „ 


This cross section is based on arguments that involve the 
conservation of angular momentum in repulsion and the assumption 
that the probability of no transition at the crossing of reactant 
and product potential surfaces increases linearly with collision 
energy in excess of the activation energy. The results are roughly 
similar to an empirical function proposed by Lotz 8 and recommended 
by Park 6 
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Eq. (7) 


In (E/E 0 ) 
s a ' E/E 0 


This function has a shape similar to experimental measurement of H 
atom ionization 7 and N atom ionization 8 by electron impact. However 
both of these functions peak at much lower energies than required 
to agree with molecular ionization cross sections by heavy particle 
impact measured by Utterback 9 and Utterback and VanZyl 10 . 

Perhaps the most complete model of reactive cross sections 
presently available is due to Landau 11 and Zener 12 . A brief review 
of their theory will establish the basis for subsequent analysis in 
this paper. Transition from a reactant system to a product system 
is presumed to occur within a narrow range of intermolecular 
distance r D where the potential surfaces of the reactant and product 
systems cross, or at least come very close to one another. The 
situation is diagrammed in Fig. (2) for repulsive potentials leading 
to endothermic reaction. The potential of the reactants at 
infinite separation is taken as the reference zero and the 
potential at the crossing point is E„. This will also be the 
activation energy E* if the potential surfaces are repulsive as 
shown on Fig. (2). However, if the potential of the product system 
should have a minimum, E* can then be greater than E 0 . The initial 
collision energy E is conserved and after the reaction the kinetic 
energy of the product species is E less the change in zero point 
energy of the chemical species involved. The product system 
potential is often only weakly repulsive, in which case the change 
in zero point energy, E°, and the activation, E*' are about equal. 

If q is the probability that transition does not occur as the 
system crosses the reaction zone, (1— q) is the probability that 
transition does occur. Then q(l-q) is the probability that 
transition does not occur during the incoming crossing but does 
occur outgoing. Similarly (1— q)q is the probability that 
transition occurs incoming and is not undone outgoing. Thus the 
total probability of transition during the collision is 

P - g(l-g) + (l-g)g - 2g(l-g) Eq.(8) 


Landau and Zener develop an exponential expression for q 


g 


exp 


( 


u fi 


~ ff 12 2 


e -a/u 


Eq. (9) 


5 


where u is the velocity with which the system crosses the reac ion 
zone and the matrix elements are the total Hamiltonian of the 
system averaged over the initial and/or final state wave func ions 


H and H 22 are just the potential surfaces of the reactant and 
product systems respectively. The gradients of these potentials 
and the interaction energy H 12 are all evaluated at r„. In general 
neither the potentials, nor their slopes, nor the interaction 
energy are known, but they are all constants which can be gathered 
in the parameter a in Eq.(9), a constant with the dimension of 
velocity. The assumptions used in the Landau-Zener treatment leave 
some doubt about the exact quantitative results in any case, but if 
the functional relations are reliable, the constant a can be 
established by comparison with measured reaction cross sections 
and/or measured rate coefficients wherever these are available. 


The following analysis considers the special case where both 
reactant and product potentials are repulsive and spherically 
symmetric with the activation energy the same as the crossing 
potential E„. If the collision is head on, that is with^ miss 
distance b equal zero, the crossing velocity n a is (2(E-E 0 ) /n) l/2 , the 
velocity of the system in center of mass coordinates when it has 
slowed from the kinetic energy E at infinite separation to the 
remaining kinetic energy at the crossing region. However, 
practically all the collisions approach with miss distance somewhat 
greater than zero. The maximum miss distance b„ which allows the 
system to reach the reaction coordinate r D is 


(K \ 2 


O ) 



Eq. (11) 


This relation obtains because of conservation of angular momentum 
and does not depend on the shape of the repulsive interaction 1 . 
However at the maximum miss distance b, the system just grazes the 
crossing surface with zero crossing velocity. At intermediate miss 
distances the crossing velocity will decrease monotonicly from u D 
at b = 0, slowly at first, and then rapidly approach zero as b 
approaches b„. The exact shape of the crossing velocity u(b) will 
depend on the shape of the repulsive potential. Whatever the exact 
shape may be, it should be possible to approximate the probability 
q(b) with a function of the type 

q{b) _ g-a/ u(b) „ e ' a/u ° [ 1- (b/b m ) d l Eq . ( 12 ) 


where d is a constant coefficient chosen to fit the true relation. 
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The reaction cross section for the collision energy E can then be 
expressed 


S(E) - 
~ 27i b 2 m 





2g(b) [l-g(i?) ] bdb 

Q-a/u 0 _ b-1 g -2 a/u„\ 


Eq. (13) 


A slow decrease in q(b) is expected near b— 0 in which case d should 
be considerably larger than unity. Then the integral of Eq.(13) is 
approximately independent of the coefficient d 


. 2 (1 -E/E) [e~ c/{E/E °~ 1)1/2 - e ~ 2c/ 

c ° 


Eq. ( 13a) 


The constant S 0 is 7 rr 0 2 times the fraction of collisions that 
approach along the potential surface leading to transition. The 
dimensionless parameter c is related to the constant velocity 
parameter a 

c . |-^J 1/2 Eg. (13b) 


The object has been merely to establish a reasonable 
functional relation between the cross section and the collision 
energy. The constants S D and c are simply chosen to agree with 
experimental data. 

Fig. (3) shows the log(S/S c ) for the Landau-Zener form of cross 
section as a function of log(E/E 0 ) for a variety of the constants 
c equal to 1, 5, 10, 20, and 30. Wherever experimental values may 
be available, the maximum in the function can be established at the 
observed value by the choice of c; then the magnitude of the cross 
section is matched by the choice of S c . In addition, Fig. (3) shows 
the more empirical functions proposed 4 '*. These more or less match 
the maximum given by the Landau-Zener model when c is a little less 
than unity, but they fall off much more rapidly for high energy 
collisions. This difference will not affect the integrations that 
provide the rate coefficient very much as long as temperature is 
small compared with E 0 /k, but can make a difference at much higher 
temperature. Also shown on Fig. (3) are the measured, cross 
sections for NO ionization by N 2 collision*. For this reaction, at 
least, a constant c about 10 is appropriate. Utterback finds some 
structure in this cross section, a very small hump near threshold, 
which no doubt indicates another reaction path. However this makes 
a negligible contribution to the rate coefficient integral and so 
is ignored here. The difference between the measurements and the 
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theoretical model is accentuated on the log plot; a better 
indication of the comparison between the two as far as the rate 
coefficient integral is concerned is a plot of S as a function of 
E/E as shown on Fig. (4). This figure shows the kind of fit to 
experiment that can be achieved with the Landau-Zener theory; in 
contrast, the empirical function ln( S/S 0 )/(S/S 0 ) is seen to peak at 
collision energies far below the observed peak and to fall off much 
too rapidly at higher energies. The same is true of the function 
of Eq.(6). While those empirical functions may approximate 
realistic cross sections for electron impact ionization or for 
dissociation by heavy particle impact, they certainly do not 
provide good cross sections for ionization of NO by N 2 impact. The 
Landau-Zener model is a far more flexible function that can better 
fit experimental data, to say nothing of the confidence provided by 
the theoretical foundation for the function. 


DISSOCIATIVE ELECTRON RECOMBINATION NO* + e -* N + 0 

The most important mechanism for ionization in air at 
moderately high temperature (the order of 5000°K or so) is the 
collision between N and 0 atoms 

N + O NO * + e Eg. (14) 


Because of the importance of this reaction, numerous measurements 
of both the endothermic and exothermic rate coefficient have been 
made. However, the scatter in the data is severe as shown by the 
exothermic rates plotted on Figure (5). Thus some uncertainty 
still exists about the exact rates for the reaction. 


The measured exothermic rate coefficients will be of most 
interest here because these have the best chance of providing the 
correct temperature dependence for the rates, as discussed above. 
The first high temperature evaluation of the exothermic rate was 
provided by Lin and Teare 13 . They found a T' 3/2 dependence of the 
exothermic rate, and chose the constant to agree with low 
temperature data of Doering and Mahan 14 . 


k 


exo 


1 . 8 x 10 21 

<Y 3 /2 


ccj 'mol -sec 


Eg. (15a) 


Subsequent experiments 15 ' 18 confirmed the low temperature result 
rather well as shown on Fig. (5) (estimates of Whitten and Poppoff 
were based on upper atmosphere airglow measurements and probably 
are not as reliable as the laboratory measurements). However, 
except for some results reported by Daiber 20 , the measured high 
temperature rates have generally been higher than Eq.(15a). 
Thompson 21 suggested that Lin and Teare's results should be 
increased by a factor of 3 , and Frohm and DeBoer 22 observed that a 
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factor of 2 increase actually gives a better fit to Lin and Teare s 
original data. Results reported by Bascomb, et al and by Stein, 
et al 24 appear to confirm the higher value. 


Since Lin and Teare 's exothermic rate was deduced from 
measurements of the endothermic process, the functional dependence 
could be in question, as shown previously. However, Dunn an 
Lordi 25 observed the electron recombination process directly m 
expanding flow from a high temperature shock tunnel. T ?} e ^ 
confirmed the same temperature dependence but with a somewhat 
larger rate 


k 


exo 


6 .7 ±2 . 3 x 10 — cc/mol -sec 

y »3/2 


Eq. (15b) 


Extrapolation of this function to lower temperature lies above most 
of the data there. Hansen 26 argued that at low temperature the 
relation should change to a T' 1/2 dependence because of the 
dependence of the equilibrium constant on temperature. Although 
this argument was based on an endothermic rate calculated for a 
constant transition probability at the potential crossing for the 
NO + and N+O systems, which is unrealistic, the same change in the 
functional trend will be found for an approximate endothermic rate 
based on the Landau Zener model of reaction. 

Two recommendations made by critical reviewers of rate data 
are shown by the dotted lines on Fig. (5). Bortner 27 reviewed the 
literature available before 1969 and recommended 

ir „ ( 1 ■ 5 ± 2 . 9 ) x 10 — cc/mol-sec Eq.(15c) 


Blottner 27 recommended a slightly different function 


1.8 x 10 19 

T 


cc/mol-sec 


Eq. (15d) 


which has been widely used in flow field calculations 3 . 

The discussion above is sufficient to indicate a considerable 
degree of uncertainty in the rate constants for Eq.(14), in part 
due to the large number of different experimental investigations of 
the reaction. In addition, an account of the effect of non- 
equilibrium on this reaction has not been derived. In the next 
section the Landau Zener model of reaction rates will be applied to 
the problem in an attempt to provide a theoretical basis for the 
temperature dependence of the rates and to include the effect of 
vibration and electronic nonequilibrium. 
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THEORETICAL MODEL FOR DISSOCIATIVE ELECTRON RECOMBINATION 

The perturbation potentials for collision of an electron with 
an NO* ion are represented by the 3-dimensional schematic of 
Fig. (6). Potentials for the N+0 repulsion and f° r the NO 
attractive well, when the free electron is at infinity, are shown 
in the U vs. r plane, where r is the intermolecular separation 
between the N and 0 atoms or between nuclear centers in the 
vibrating NO* ion. The repulsive N+0 potential is presumed to cross 
at the ion's potential minimum; then the activation energy for 
endothermic reaction is just the change in zero point chemical 
energy and the activation energy for exothermic reaction vanishes, 
which seems to be in accord with experimental observation. The 
electron approaching the ion with kinetic energy E. is accelerated 
as it falls into the strong Coulomb attractive potential (E 0 -e 2 /r e ), 
which is drawn in the U vs. r e plane; r e is the distance between the 
electron and the ion. Eventually the electron experiences a strong 
repulsion as it encroaches on the space occupied by the ion's 
electrons. This slows the electron to its initial kinetic energy 
and opens the possibility of transition to the N+0 system if the 
vibrational configuration of the ion is nearly at the potential 
crossing point r D . Of course, electrons behave more like guantum 
waves than classical particles in collision. However, except for 
misrepresenting features like glory and rainbow effects, a 
classical treatment is expected to cut through the average of the 
quantum results. 

Because of the large Coulomb attraction, the cross section 
which leads to a close enough encounter to permit transition is 
very large. The maximum miss distance b„ which allows the electron 
to come within some distance r* needed (the order of a Bohr radius) 
varies about inversely with the original electron energy. This is 
given by Eq.(ll) except that in this case the potential is 
attractive 


JL _ (M 2 . 1 + _ef_ « fi- Eq.(lla) 

S* U*J r* E e r* E e 


The electron velocities will be the order of 100 times larger 
than heavy particle velocities in the gas, so for most collisions 
the transition probabilities are expected to be small. Then the 
total probability of transition given by Eq.(8) is approximately 
inversely proportional to the electron's velocity 


P - 2(1 - e~ a/Ua ) 


2a 

u e 


const : 

El ' 2 


Eg. ( 16) 
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Although this approximation is violated at very low collision 
energies, these are relatively unimportant for exothermic reaction 
as they are weighted so lightly by the Boltzmann distribution of 
energies. Eqs. (11a) and (16) combine to give the functional 
relation for total cross section which is inversely proportional to 
the 3/2 power of the electron's kinetic energy 


S « — P 

S' 


const : 

El ' 2 


const 

(kT e ) 3/2 x 3/2 


Eg. (17) 


Then the rate coefficient for transition has the form 


- »/ 


S x e~ x dx <* 


— f 

T„ J 


cbc 


.1/2 


const 

T„ 


Eg. ( 18) 


Eg. (18) gives the functional dependence of the rate 
coefficient for a favorable vibrational configuration of the NO* 
ion. The probability p v that an ion in vibrational state v is 
within a small region <5 near the bottom of its potential well, and 
thus within the reaction zone, is 


6 _ / n \ 1/2 6 

T u v \2W5 to) T 


Eg. (19) 


where r is 2ir/ts>, the period of vibration. The overall rate 
coefficient is obtained by summing the rate coefficients for each 
state v weighted by the probability that the ion is in the v* K 


vibrational state, as in Eg. (5b) 


k « 


1 

O T 

v e 


E 


(v 0) 1/2 


! (nT v ) 1/2 
0 v T e 0 


Eg . ( 20a ) 


Then approximating Q v with the harmonic oscillator partition 
function 

k « ( i - e~ 0/Tv ) Eg. (20b) 

T 


k 


r v >e > 


const 

rrt rpl/ 2 

1 e 1 v 


T =T =T 

x e v 


cons t 

fj 1 3/2 


Eg. (20c) 
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At very high temperature the overall rate coefficient is 
predicted to vary inversely with the electron temperature and 
inversely with the square root of the vibrational temperature. In 
many practical cases T e and T v will be nearly equal, and if fu 
equilibrium obtains, the overall theoretical rate will vary as T 
in accord with the experiments of Dunn and Lordi 25 . The constant is 
chosen so that at high temperature equilibrium the value will match 
Dunn and Lordi's measurement, Eq.(15b); then the rate coefficient 
for dissociative electron recombination with NO* is 

1/2 

k - 3.0 X 10 18 — — ( 1 - e~ 3420/7V ) cc/mol- sec Eq.(21) 


The temperatures are in °K. At very high temperature the 
vibrational partition function and the sum over vibrational states 
in Eq.(20a) should be truncated and anharmonic effects included. 
However, since these appear in ratio with one another these 
corrections will tend to be compensating. Since this reaction is 
exothermic, the effect of excited electronic species with different 
vibrational frequencies can be discounted, since the small 
populations in these excited states have little effect on the total 
rate observed. Thus Eq.(21) is assumed to be a relatively reliable 
rate coefficient for the dissociative electron recombination with 
the NO* ion. This function will next be used with the reaction's 
equilibrium constant to determine the endothermic rate coefficient 
for ionization resulting from collision of N and O atoms. 


IONIZATION DUE TO N + O COLLISION 

If the exothermic rate coefficient for NO* and electron 
recombination is known, the endothermic rate for ionization due to 
collisions between N and O atoms can be obtained by dividing the 
exothermic rate with the equilibrium constant. However, in the 
present case, a state of pseudo-equilibrium is considered where the 
electron temperature, the vibrational temperature, and the heavy 
particle kinetic temperature are approximately decoupled and 
frozen. Then the exothermic rate is presumed given by Eq.(21), 
while the equilibrium constant needs to be evaluated for the 
pseudo-steady condition rather than at full equilibrium. 

A pseudo-steady state of a chemically reacting mixture 
requires that the change in chemical potential created by the 
reaction must nearly vanish 

Ap - H*,* + He + + * 0 

Eq . ( 22a ) 

- A E° -kT In Q(NO + ) -kT 0 In Q(e) 

+ kT In Q(N) + kT In 0(0) 


12 


Eq. (22b) 


where the change in zero point energy is just 


A E° 


NO' 


- E‘ 


N 


- E 


O 

O 


since the zero point energy of the free electrons is taken as zero. 
Then 

A 2? ° , Q(NO') [C>(e)] Tjr Eq.(23) 

kT " Q(N) 0(0) 


This pseudo- steady relation differs from the full equilibrium 
relation in that the electron partition function is raised to the 
power T e /T due to the fact that the electrons have a different 
kinetic temperature than the heavy particles. 

The partition functions can be either pressure standardized or 
concentration standardized 29 . Since the number of mols is constant 
in this reaction, the equilibrium constant is the same in either 
case. The notation Q p (i) and Q c (i) will be used for the pressure 
standardized and concentration standardized partition functions for 
specie i, respectively. 


A E° 

kT 


Q p (NO+) [Q p (e)] TjT 
Q p m Q p (0) 

Q c (NO*) [Q c (e) ] TJT 
11 Q C (N) OJO) 


In 


In 


P W o> tPj Ta/r 
Pn Po 

n N( y [nj TJT 


n N n o 


Eq. (24) 


where and n t are the partial pressures and densities of species 
i in the chosen standard units (typically atmospheres pressure and 
either molecules/cc or mol/cc, respectively). In concentration 
units, the equilibrium constant can be expressed 


n N0 ' n e q . Q c ( NO *} 0c ^-\E°/kT Eg. (25) 

n N n Q ’ Q C (N) Q c (0) 


This has the same functional form as the full equilibrium constant 
except for the correction factor G 


G 


Oje) j 


Eq. (26) 
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which is a function of the electron concentration n e and the 
temperature ratio T./T. G becomes unity when the kinetic 
temperatures are equal, of course. 

To the first order approximation the upper electronic states 
are neglected and the diatomic species are considered rigid 
rotating harmonic oscillators. Then the functional relations of 
the concentration standardized partition functions are in standard 
units of mol/cc 

Q(NO + ) « ^—T7t ~ Eq . ( 27a ) 

(l-e~ e/7V ) 


(?(e) « tI /2 

Q(N) « T 3/2 

Q(O) « T 3/2 


Eq. (27b) 
Eq. (27c) 
Eq. ( 27d) 


To this order, the functional expression for the equilibrium 
constant is 


- 3/2 


*c - 


endo 


-A E/kT 


T 1/2 (l-e“ 6/TV ) 


Eq. (28) 


and the form for the endothermic rate coefficient becomes 


^ endo 


k K <* G 

exoc 


rjiXi 2 fTil/ 2 

1 V 1 e___ e ~bE/kT 


T v -T e -T > 0 


^-. 1/2 

> const t 1 ' 2 e~* B/kT 


Eq. ( 29) 


The high temperature equilibrium limit is the same relation 
proposed by Blottner 28 and used by Gupta, et al 3 . 

The constant of Eq.(29) is determined using quantitative 
values of the partition functions. Again to first order, 

neglecting upper electronic states and considering the diatoms as 
rigid rotator harmonic oscillators, the logarithm of the 

concentration standardized partition function is for the molecules 

ln£> -In ( 7 + In f-X) - In (l-e‘ e ’' /r '') + In g 0 

\ h 2 ) V®*/ Eq.( 30a ) 

- — In T + — In M -8 . O725-ln0 r -ln (l-e' e, ' /Tv ) +ln g Q 
2 2 
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Eq. ( 30a) 


while for the atoms it is 

In O - —In T + —In M - 8.0725 + In g Q 
v ° 2 2 


where M is the molecular weight in atomic mass units, 0 r the 
rotational energy constant B/k in degK, 0 V the characteristic 
vibrational temperature ho/k in degK, and g 0 the degeneracy of the 
ground state, or low lying states, of electronic excitation. For 
the species of interest here the concentration standardized 
partitions functions in units of mol/cc are 


In Q c (NO+) 


— In T - 4 . 0286 
2 


In (1 - e' 3 * 2a/Tv ) Eq. ( 31a) 


In 


— In r. - 19. 3324 
2 e 


Eq. (31b) 


In Q.AN) - — In T - 2.7276 


Eq. (31c) 


In QAO) 


— In T - 1 .7164 


Eq. (31d) 


and the change in zero point energy is 

A E°/k - 32400 °K 


Eq. ( 31e ) 


Accordingly, the first order approximation for the equilibrium 
constant is 


,3/2 


K - 1.218 x 10~ 8 G • 


t i /2 {1 _ e -3420/r v) 


e -3240Q/r Eq.(32) 


and the endothermic reaction rate which results from multiplying 
Eq .(28) with Eq.(21) is 


k 


endo 


4 . 8 x 10 10 G ■ 


tjjXI 2 
1 v 1 e 

rpl/2 


e -32400/r 


cc/ mol -sec 


Eq. (33) 
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At high temperature equilibrium this rate approaches 


^ endo 


T s -T v -T>d 


> 4.8 x 10 10 T 1/2 e~ 32i00/T cc/mol-sec Eq.(34) 


For concentration units in mol/cc, the free electron partition 
function of Eq.(31b) is used and the correction factor G may be 

expressed 

log G - (1 -TjT) ( log n e - -^-log T e + 8.396) Eq.(35) 


Fig. (7) shows the ratio of log G to (1-T e /T) as a function of 
electron temperature T e for electron concentrations of 10 , 10 , 

10' 9 and 10' 11 mol/cc. The highest density corresponds to about 10-s 
ionization at a total density of one amagat, and the smaller 
densities correspond to smaller degrees of ionization or to smaller 
total density at higher altitude, for example. 


HIGHER ORDER PARTITION FUNCTIONS 

At very high temperature some additional terms in the 
partition functions become important because the upper electronic 
states are excited and the effects of vibration-rotation coupling 
and the anharmonicity of the diatomic species become more 
significant. Jaffe 30 shows that more rigorous partition functions 
for these species are duplicated reasonably well by the Dunham 
potential approximation, for which the rotational-vibrational 
energy of the diatomic particle is expressed in terms of the 
quantum numbers j and v respectively 

- BjU* 1 ) -(^ 

+ u>v - (<jix) v 2 - a vj ( j + 1 ) 


j 2 (j + 1) 


Eg . ( 36 ) 


Note that the older spectroscopic notation is used in Eq.(36), 
which is sufficiently accurate for the purposes at hand, and that 
the second rotational energy constant D (no relation to 
dissociation energy) has been given the value (4B 3 /w 2 ) as derived by 
Pauling and Wilson 31 and Herzberg 32 . The expression for the 
partition function can be expressed 

0 - 0 C £ 9t e~ ei/kT Q r j 0 vi d + r,T) Eg. (37) 


where the factors Q t , Q ri , and Q vl are the usual separable partition 
functions for translation, rotation, and vibration respectively. 
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The summation is over all electronic states i with energy e t and 
degeneracy g, ; Q„ and Q„ come under the summation since the 
rotational and vibrational constants all depend on the 
state, but Q t does not. 

Mayer and Mayer 33 derive the factor r which corrects the re sult 
for vibration-rotation coupling and anharmonicity for the case 
where the temperatures are all m equilibrium. In the present 
problem the electronic and vibrational temperatures are taken eq 
to T and rotations are considered to be at equilibrium with the 
kinetic gas temperature T. The correction factor rederived for 
this state of nonequilibrium becomes 


1 + T 


1 + 


*BA kT + f «i 

J V B i 


2 (o)x) kTy Eq . (38a) 


w i 




The rotational and vibrational constants for the different 
electronic states have been tabulated in convenient sb< p rt f °*' m ¥ 
Jaffe 30 and Park 6 . One additional correction that should be included 
at high temperature is truncation of the rotational and vibrational 
partition functions at the dissociation energy D so they do not 
become infinite as they would for and infinity of states. The 
separable rotational partition function then is 


kT 

B 4 


( l - e' Dl/kT ) 


Eq. (38b) 


while the separable vibrational partition function is 


Eq. ( 38c) 


The ratio of the improved partition function for N0+, Eq.(37), 
to the first order expression Eq.(30a) is graphed in Fig. (8) for 
three ratios of vibration-electron temperature to kinetic 
temperature: 1.0, 0.6, and 0.3. At 50,000°K the more ex ^ct 

partition function can become 5 to 10 times larger than the 
simplified low temperature expression, though diatomic species may 
not last long at that temperature, of course. 

The large multiplicity of electronic states of the atoms near 
the ionization limit poses a problem; the partition functions grow 
to infinity if all the states are included, but the effective 
cutoff for these states is not so clear. Often the ionization 
limit is lowered 34 to the level where the free space available o 
the outer electron is limited by the density of the gas, °r b y the 
Debye shielding effect of electrons surrounding the ion which 
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truncates the range of the ion's Coulomb potential, or by 
collisions with electrons which broaden the energy levels until 
they effectively become part of the continuum for the free 
electron. These criteria are awkward to use because they depend on 
different properties of the gas mixture 29 ; also the cutoffs are 
generally at such high guantum numbers that the partition functions 
are still suspiciously large. The table of electronic energy 

levels for the atoms provided by Park 6 terminates at quantum level 
10, which might be sufficiently close to the ionization limit for 
these almost exactly hydrogen eigenvalues. A simple cutoff will be 
used here which should apply to the case where the electron density 
is far below the equilibrium value at least; the upper electronic 
states within about kT of the ionization limit I will be depleted 
by rapid escape to the ionized continuum with few reverse reactions 
to replenish these upper states, leaving a pseudo steady decreased 
population in these states 351 . In order to provide a smooth 

transition of the upper state populations, the equilibrium 
distribution is multiplied by an effective cutoff function 

Jk = (i - e - {I - ei)/kT ° ) gi e~ ei/kT * Eq . ( 39 ) 

n 

Fig. (8) also shows the calculated ratios of the improved 
partition functions to the first order functions for the N and 0 
atoms, again for three different ratios T e /T. 

The final correction to the equilibrium constant is the ratio 
of Q ( NO* ) /Q 0 ( NO* ) to the product of Q(N)/Q 0 (N) and Q(0)/Q o (0) from 
Fig. (8). This correction factor is shown on Fig. (9). Due to the 
compensation provided by these ratios, the . first order 
approximation for the equilibrium constant is within a factor of 
two of the more exact values at full equilibrium; the correction 
rapidly approaches unity as the vibration electron temperature 
decreases. In view of the approximations inherent in the model, a 
further correction of this small amount seems unwarranted. 


HEAVY MOLECULE COLLISION IONIZATION OF NO 

The ionization of NO by collision with heavy molecules in the 
gas is so infrequent that this reaction mechanism is usually much 
less important at high temperature than the collision of N and 0 
atoms, although this reaction is sometimes included in schemes for 
solving flow of chemically reacting air 3,28 . However, this reaction 
is particularly interesting in view of the fact that it is one of 
the rare cases where measured collision cross sections 9,10 can be 
compared with an accepted rate coefficient. 

The cross section function can be expressed 

— - qr(l-qr) Eq.(40) 

5 o x 
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where q is the Landau Zener probability that no transition takes 
place at the potential crossing region 


g(y) = e 


-c(x„/y ) 1/2 


Eg. (41) 


The rate coefficient of Eq.(l) then becomes 

kf - u S 0 F(x 0 ) e~ x ° 


Eq. (42) 


where the Arrhenius preexponential factor F(x„) is 


F(x 0 ) - 2 J q (1 -q) y e' y dy 

O 


Eq . ( 4 3 ) 


This inteqral was performed by numerical quadrature and the results 
Ire given 9 ” Fig. (10) where the logarithm of F(x„) is shown as a 
function of the logarithm of (c’x.). For the NO ionization cross 
sections shown on Fig. (4), the constant c is about 10. When « 
reasonably sized total cross section S 0 of 10 cm is used, t 
resulting rate coefficient is much smaller than the accepted rate 

coefficient 3 - 28 


k f “ 


2.2 X 10 15 e -108000/T 


i.35 


cc/mol-sec 


Eq. (44) 


This occurs because the excited states are not included in the 
computation, as they should be at high temperature . The total ra 
is a summation over all excited states as in Eq.(5b) 


*f - ' 


uS„ 1 

O 


^-\ f* 2 (Xj) 

V (l-e i f I ) 2 


Eq. (45) 


The cross sections are increased by (1-ei/I)' 2 , where e t is the 
electronic energy of state i, to account for the increased size of 
the excited state wave function, which in the upper states at leas 
is almost hydrogen like. The summation is approximated for all the 
electronic, vibrational, and rotational states available to the NO 
molecule. The resulting ionization rate coefficient is shown by 
the dashed curve of Fig. (11), where the reference total c f os s 

_ . w-w -is 2 *i c? ar«h i D\ 70 n Pi I 


section S c is chosen 1.6 x 


10 


cm' 


so agreement is achieved at 


10000°K. Note that this factor includes the fraction of collisions 
which approach along the transition potential. The principal point 
is that the functional dependence predicted by Eq.(45) can be 
brought into reasonable accord with the simple Arrhenius expression 
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of Eq . ( 4 4 ) . 
extrapolation 
temperatures . 


This lends some extra justification f° r the 
of the Arrhenius function to much higher 


CONCLUDING REMARKS 

The Landau Zener theory of reaction cross section has been 
used to derive a functional expression for dissociative electron 
recombination with the N0 + ion. The rate coefficient is predicted 
to vary inversely with the electron temperature and the vibrational 
partition function and directly with the square root of vibrational 
temperature. At full equilibrium these temperatures equal the 
kinetic temperature, in which case the rate coefficient varies 
inversely with the square root of temperature at normal 
temperatures and inversely with the 3/2 power of temperature at 
temperatures that are large compared with the characteristic 
vibrational temperature, -fiu/k. This relation is more or less in 
agreement with most of the experimental determinations made at both 
normal and high temperatures, whereas the simple T' V2 variation that 
has usually been used when making flow field calculations is unable 
to simultaneously describe closely the observations over the entire 
range of temperature. Since many of the flow fields around the new 
generation of space vehicles will involve highly nonequi librium 
vibrational and electronic distributions, the specific variation of 
the rate coefficient with each of these temperatures is required. 
The present model is able to address this need in an approximate 
manner, at least. 

Unfortunately, measurements of cross section are not available 
for dissociative electron recombination with ions. Thus, 
comparison with the observed rate coefficient is the only way that 
the appropriateness of the Landau Zener theory to these reactions 
can be judged. However, in the case of a reaction of relatively 
minor importance, heavy particle impact ionization of the NO 
molecule, measured cross sections are available, at least for the 
ground or low level vibrational and electronic states of the 
molecule. The Landau Zener model can be adjusted to fit these 
measured cross sections reasonably well and then integrated over a 
Maxwell Boltzmann distribution of collision energies. The 
resulting rate coefficient applies only to ground state NO and so 
is much smaller than the total rate coefficient for a gas which 
includes higher states of the molecule. However, when the excited 
states are assumed to have the same functional form of cross 
section and the total scattering size is taken as the size of the 
outer electron wave function, the summation over all excited states 
provides a rate coefficient in reasonable accord with the equation 
that has been used to represent the total rate coefficient for 
heavy particle collision ionization of NO. 
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0 2 DISSOCIATION BY Ar COLLISIONS 


SCHEMATIC DIAGRAM OF REACTION POTENTIAL SURFACES 


REACTION CROSS SECTIONS S/S 0 AS A FUNCTION OF 
COLLISION ENERGY E/E 0 


REACTION CROSS SECTION FOR IONIZATION OF GROUND STATE 
NO BY N, COLLISIONS 


EXPERIMENTAL AND THEORETICAL RATE COEFFICIENTS FOR 
DISSOCIATIVE ELECTRON RECOMBINATION WITH NO IONS 


SCHEMATIC DIAGRAM OF POTENTIAL SURFACES FOR 
e, NO* , AND N+O SYSTEMS 


Fig. (7) 


Fig. (8) 


CORRECTION FACTOR FOR EQUILIBRIUM CONSTANT EXPRESSION 
WHeS APPLIED TO PSEUDO-STEADY NONEQUILIBRIUM ELECTRON 
TEMPERATURE CONDITIONS 


RATIO OF CORRECTED PARTITION FUNCTIONS TO FIRST ORDER 
EXPRESSIONS FOR NO*, N, AND O 


Fig. (9) 

Fig. (10) 
Fig. (11) 


RATIO OF CORRECTED TO FIRST ORDER EXPRESSION FOR THE 
EQUILIBRIUM CONSTANT OF THE N + O - NO* + e REACTION 


PREEXPONENTIAL FACTOR FOR LANDAU ZENER RATE COEFFICIENTS 


ENDOTHERMIC RATE COEFFICIENT FOR NO + N 2 -* NO* + e + N 
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Fig. ( 1 ) 


Oj DISSOCIATION BY Ar COLLISIONS 
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Fig. (4) REACTION CROSS SECTION FOR IONIZATION OF GROUND STATE 
NO BY N 2 COLLISIONS 
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Fig. (11) ENDOTHERMIC RATE COEFFICIENT FOR NO + N 2 


NO + + e + N 2 


